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Abstract 

We enumerate weighted graphs with a certain upper bound condition. 
We also compute the generating function of the numbers of these graphs, 
and prove that it is a rational function. In particular, we show that if 
the given graph is a bipartite graph, then its generating function is of the 
form rj~TSTi > where m is the number of vertices of the graph and p(x) 
is a polynomial of degree at most m. 



Key words Weighted Graphs, Rational Convex Polytopes, Rational Gen- 
erating Functions, Ehrhart (Quasi-)polynomial. 

1 Introduction 

For a given nonnegative integer n, let [n] := {0, 1, 2, • • • , n}. Also, let G = (V, E) 
be a simple graph (no loops and no multiple edges allowed) with the vertex set 
V = {vi, v 2 , • ■ ■ , v m }. Let a = (ni,n 2 , • ■ • , n m ) £ [n] m so that 



•j 



< n if ViVj G E. (1) 



In other words, the sum of two weights corresponding to an adjacent pair of 
vertices is bounded by a given integer n. 

Wc call a triplet WG a — (V, E, a) a weighted graph of G with distribution 
a. We will denote by WG(n) the number of all weighted graphs of G with a 
fixed upper bound n. Note that WG(n) enumerates the set of all solutions of 
the system of linear inequalities corresponding to the graph G by the condition 
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CD- 

Let Q be the set of all simple graphs and C[[x]] the ring of formal power 
series. Define a map 

oo 

p : G -► C[[x]] (G i ► p(G) = WG(n)a; n ). 
For example, consider the following graph G: a c . 



d e 

For each a = i G [n], b,d,e can be 0, 1, 2, • • • , n — i, and for each b = j G 
[n — i], c can be 0, 1, 2, • • • , n — j. Hence, 

n n~i 

WG(n) = J^in + 1 - *) 2 (E( n + 1 - •?)) 

i=0 j=0 

_ (3n + 5)(2n + 3)(3n + 4)(n + 2)(n + 1) 
~~ 120 



and 



1 . „ 1 + 26a; + 66a; 2 + 26a; 3 + x 4 „ 1 + 11a; + 11a; 2 + a; 3 

P{G) = T20 (18 + 45 O^F 

+ 40 1 + 4 * + f 2 +15 l + X 



(l-x) 4 ' ^(1-x) 3 ' 2 (l-a;) 2 ^ 
(1 + x)(l + 7x + x 2 ) l + 8x + 8x 2 + x 3 



(l-x) 



(l-x) 



Equivalcntly, WG(n) is the number of solutions (n a , rib, n c , rid, n e ) G [n] 5 to 
the following system of inequalities: 

n a + nt, < n 
rib + n c < n 
n a + n d <n 
n a + n e < n 



Remark 1.1. 

1. The map p is not injective. 



2 



For graphs G\ * » and Gi 



9 n+l n+l 

l + Ax + x 2 



fe=i fc=i 

PFe -used i/ie following identities (to get p{G) in the previous example): 



(n + l) m = Mm, k)C(n + k, 



k=l 



where A(m, k) is an Eulerian number, and C(n, m) = m \^- m )\ * s a bino- 
mial coefficient (See Fact \2.5i 3. in Section WIK ). and 



oo ^ 

J2C(n + k,k)x n = {l _ x)k+1 



n=0 



In Section [2] we compute p{G) for linear graphs, circular graphs, complete 
graphs, star graphs, discrete graphs, a cubic graph, an octahedral graph and 
complete bipartite graphs. In Section [3] we remind a reader of rational con- 
vex polytopes and rational generating functions, and describe the relationship 
between these objects and our problem (enumeration of weighted graphs). In 
Section 0] we list many new problems that we encountered while we worked on 
this topic. 

2 Generating Functions for Various Graphs 
2.1 Linear Graphs 

For the null graph we set W$(n) = 1 for all nonnegative integers n, so 

1 



1 - x 



Let L\ be a one-vertex graph, and for i > 2 let Li be the linear graph with % 
vertices. That is, i,; is a tree with two vertices of degree one and i — 2 vertices 
of degree of two. For the one-vertex graph L\ = Ki(»), WLi(n) = n+l and 

p{Li) = jihef- 

For L2 = K2 ( a» # ^ ) let a = i, b = j, and 
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, I 1 if i + j < n, . 2n + 3 — i — j . 

Oij = < = I J. 

I otherwise n + 2 

Note that all indices i and j run from to n, not from 1 to n + 1. B(n) := (by) 
is an (n + 1) X (n + 1) matrix, J„ := (1, 1, 1, • • • , 1)* e [n] n+1 . Then 

n 

WL 2 (n) = JlB{n)J n = ^ b tj = C(n + 2, 2), 

i,3=0 



For i 3 ( . . . ), 



WL 3 (n) - J*B(n) 2 J n = 
P(L 3 ) = 

By induction on k we can show that 



C(n + 3,3) + C(n + 2,3), 
1 + x 



(1 



holds for all nonnegative integers /c. 
For example, 



B(4) 



1111 
1111 



1110 
110 
1 



B(AY 



1111 



15 


14 


12 


9 


• r j 


14 


13 


11 


8 


4 


12 


11 


9 


6 


3 


9 


8 


6 


4 


2 


5 


4 


3 


2 


1 



B(A) 4 



Note that WL k+l {n) is the (0, 0)-entry of B(n) k+2 . For example, WX 3 (4) 
55, which is the sum of all entries of the matrix B(4) 2 . 



55 


50 


41 


29 


15 


50 


4G 


38 


27 


14 


41 


38 


32 


23 


12 


29 


27 


23 


17 


9 


15 


14 


12 


9 


5 
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Theorem 2.1 (Bona and Ju, 2006 [8j). Let 

oo oo 

F(n,x) = Y,WL k+1 {n)x k = ^(J*S(n) fc J„; 

k=0 k=0 

Then 

1 + F{n - 1, -x) 1 



F(n,c 



■r 



l-x(F(n-l,-x) + l) 1 



1 + F(n - 1, -x) 
-x, 1, a;, 1, F(n — 2, x)], 



l + F(n-2,x) 

where F(0, x) = ^ = [-x, 1] and F(l, 0) = = [_ a!j 1>a?i i]. 

Theorem 2.2 (Bona and Ju, 2006 |5|). 

7-1/ \ Pn\Xj 

= 04^)' 

Q n (x)=det(/-xB(n))=^ L 2 J (-l)^J^ 

fc=o ^ ' 

and 

Pn{x) = (Qn-2{x) - (1 + X)Q„(X))/X 2 

is a polynomial of degree deg{Q n (x)) — 1 = n. 



Remark 2.3. 

1. Q n (x) satisfies a recurrence relation 

Q n {x) + {x 2 - 2)Q n _ 2 (x) + Q n -i{x) = 0. (2) 

2. From the recurrence relation (J2J we get the following generating function 
R(x,y) for the sequence {Qn(x)}^L Q : 

T,, x V^-n / \ n {I + y){l - xy - y 2 ) 
^ 1 -y 2 {2- x 2 -y 2 ) 
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3. {F (n, X )}^ =0 is an approximation(called Fade approximant) to the infinite 
periodic ( of period 4 ) continued fraction 



X X ^ IT/ > X • X ■ X • »37 ^ X • — 



and converges to 



2x 

-(2 + .)W^4 ^ [2)0o)j 

on (-00,2]. 



^. The explicit expression of the generating function G{x,y) for the sequence 
{F(n,x)} < ^L is not known, but from Theorem 2.2 



oo oo „ / ^ 



1 + x 



lL -' n x 2 Q n (x) x 2 (l-y) 

n—O n— v ' v 7 

5. Lei s(M) 6e i/ie sitm o/ all entries of the mxm matrix M with nonnegative 
integral entries, and let M^ k \x) be the mxm matrix I — xM where the 
k-th row of the matrix I — xM is replaced by (1, 1, • • • ,1). fk(x) be the 
determinant of the matrix M^ k '(x). Then the generating function r](x) for 
the sequence {s(M n )}^L is as follows (the proof is immediate from the 
4.7.2 Theorem of R. Stanley JES) ■' 



oo \r-^ rn r / \ 

V (x) = Y s(M n )x n = }^ l = , L 
det(I-xM) 



2.2 Circular Graphs 

Before considering circular graphs, we first consider the linear graphs Lk+\ with 
k + 1 vertices, two of them being end vertices a and 6, as follows below: 
a* • • ■ • ■ • • • u b 



If we fix n a — i and nb = j, then the number of all possible ways to give 
distributions corresponding to the vertex set of the graph Lk+i is same as the 
entry of the matrix B(n) k . 

Next, we identify the leftmost vertex a and the rightmost vertex b, resulting 
in the circular graph Ck with k vertices. 

Hence, the number of all possible ways to give distributions corresponding 
to the vertex set of the graph Ck (with identified vertex having n a = = i 
fixed) is same as the (i,i)-entry of the matrix B(n) k . 

This implies that WC k {n) = tr{B(n) k ) for k > 2, and Wd{n) = trace(I k+1 ) = 
n+1. 
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Lemma 2.4 (Theorem 4.7.2 and Corollary 4.7.3, Stanley, 1997 |2T]). 

For any m x m matrix M , we have: 

^ x-f (det(7-xM)) 

^ v ; det(/-xM 



By the Lemma \'2. 41 above 



CF(n,x) = f>(B(»)V = \?±±\z- 

and 

We list several p(Ci)'s and ,o(L;)'s for reference. 

P(Ci) = 7i=SJ» = ti=SJ» 

P(C 2 )- IT ^ F ^ p(i 2 ) = IT ^ F 

p(Ces) = ( i+tf(|f a; )4 p{Lz) = (jzfy? 

Pl°4j - (l_ x )5 p[pi) - (l_ x )5 

„/r< \ _ l+6a: + llx 2 +6a 3 +a: 4 r \ _ l+7x+7a; 2 +a: 3 

^°5-> - — (i+x)(i-x)* — Py^) - — (T^F — 



For example, WC 3 (n) = jg (4?i 3 + 18n 2 + 28n + 15 + (-1)") 
f |(n + 2)(2n 2 + 5n + 4) if n is even 
l)(2n 2 + 7n + 7) if n is odd ' 



2.3 Discrete Graphs 

In order to describe discrete graphs, complete graphs, star graphs, complete 
bipartite graphs, etc, we need the notion of Eulcrian numbers and Eulcrian 
polynomials. 

Let p = P1P2 ■ ■ ■ Pt be a t-permutation. We say that i is a descent of p if pi > 
Pi+i - Let A(t, k) be the number of t-permutations with k — 1 descents. The num- 
bers A(t, k) are called the Eulerian numbers, and A t (x) = J3k=i ^(^' ^) a;fe is 
called the Eulerian polynomial. 



Fact 2.5. FFe list below several known facts about Eulerian numbers and Eu- 
lerian polynomials. [See Bona, 2004 an d 2005 ^ for details about Eulerian 
numbers, and also Graham et. ai, 1994 03/ 

1. A(t, k) = kA(t -l,k) + (t-k + l)A(t -l,k- 1). 

2. A{t, k) = A(t, t+l-k) and YX=\ A (t, k ) = t-- 
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3 - x * = ELi A (t, fyC{x -1 + k, t), where C{t, k) = fc!( /l fc)! . 

I g{x,y) = EZo ELi A (t, k)x% = TZoM^i = iz^t - 

c A t (x) _ d r A t _i(x) 1 
^' xil-xy+i ~ dxl (l-x)* J" 

6. A t {x) = teA t _i(x) + x(l - x)A' t _ 1 (x). 

7. A t (x) = x(l ~ xf +J2l =1 C(t,k)(l - xy~ k A k (x). 

8. A{t,k)=Yll= {~V l C(t + l,i){k~i) n - 

Let D t = (V, E) be a discrete graph, that is, a graph with 1^1 = t and E 
Then, clearly, 



WD t (n) = {n + l)\ P {D t ) 



{1-XY+ 1 ' 

and 



F(x l y)=X;p(A)^ = r 



2.4 Complete Graphs 

Let K t be a complete graph of order t. Then, 

Theorem 2.6 (Ju, 2006 [T6]L 

r=l 

and 

OO OO £ y 

F(x, y) := ^ E ^W^fr = (! + *+ TZ^T^r 



We define 



a(y; x) = 1 + x + - , (3(y,x) 



1 — x ' ' 1 — a; 2 e !/ 

Using the facts given in Section 1231 we can show that the following holds: 



P[ y, x ) - 1 _ x2(iy - ^ t! )( 1 _ x2 + ^ (1 _ ^ 2)t+1 t] j 

i-x 2 ^ x 2 (i-x 2 y+ 1 1\ ' 
*=i v y 
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and 

F(xy) - 1 , V -y A t {x 2 )+txA t -i{x 2 ) y t 
{ ,V) l-x + (1-x) 2 + ^ x 2 (l-x)(l-a; 2 ) t t! * 1 j 

For t>2, the numerator of the summand in the equation J^J is 

Atix^ + txAt^x 2 ) .txA t -i(x 2 )+x 2 (l-x 2 )A' t _- L (x 2 ) 
-2 = ( 1 + x ) ~2 ■ 



Let 



. . txA t _ 1 (x 2 )+x 2 (l-x 2 )A / t _ 1 (x 2 ) 

rt ( X > = ~2 • 



for t = 2, 3, • • • . By Fact 12.51 2. (symmetric condition), 



t-i 

r t (-l) = Y,( 2k ~ *)^(* - 1> fc ) = °' 



fe=i 

for t = 2, 3, • • • . This implies that 

A t (x 2 ) + txA t „ 1 {x 2 ) 



(1 + x) 2 p 2t -i{x) 



for some polynomial p2t-i{x) in x of degree 2£ — 4. 

Theorem 2.7. for C/ t (x) = J2n=o WK t {n)x n , 

T7T/ \ tt I \V 1 V , ST V2t-i{x) y 

where 



x (l_ x )2 (1 _ x) 3 (1 _ a .2 )t -2 t ,. 



P2t - 4(a;) = ^ 2 (TT^) 2 

is a polynomial of degree 2t — 4 /or i = 2, 3, 4, • • • . 
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Computations (using Maple) of Ut(x) for t = 0, 1, 2, • • • ,8 have shown that: 
1 

l 

(1 -x) 2 ' 

1 

(1 - x) 3 ' 

1 + .x + x 2 
(l-x) 3 (l-x 2 )' 
1 + 2x + 6x 2 + 2x 3 + x 4 

(l-x) 3 (l-x 2 ) 2 
1 + 3x + 19x 2 + 14x 3 + 19x 4 + 3x 5 + x 6 

(l-x) 3 (l-x 2 ) 3 
1 + 4x + 48x 2 + 56x 3 + 142x 4 + 56x 5 + 48x 6 + 4x 7 + x 8 

(1 -x) 3 (l -x 2 ) 4 
1 + 5x + 109x 2 + 176x 3 + 730x 4 + 478x 5 + 730x 6 
(l-x) 3 (l-x 2 ) 5 
176x 7 + 109x 8 + 5x 9 + x 10 
+ (l-x) 3 (l-x 2 ) 5 

Remark 2.8. 

1. The sequence a\,a2,- ■ ■ , a n of positive real numbers is called unimodal if 
there exists an index k such that 1 < k < n, and a\ < 02 < ■ • • < afc > 
flfc+i > • ■ ■ > o-n- The same sequence is called log-concave if ak-iau+i < 
a 2 holds for all indices k. It is well-known that a log-concave sequence is 
unimodal. (See Bona\ljj for the proof.) The sequence {A(t, k} f k=1 of Eule- 
rian numbers is log-concave, so it is unimodal for all t. The coefficients 
of numerator in Ut(x) is unimodal for t = 1,2,3,4,6, but not always as 
shown in C/s(x), Ur(x) above. 

2. The rational functions Ut(x) have denominator (1 — x 2 ). This means 
that WKt(n) is a (in fact, Ehrhart) quasi-polynomial of a certain poly- 
tope. (We will mention this later in Section^) Hence, its form depends 
on the parity of n. (See the next remark.) 



U {x 
Ui{x 
U 2 (x 
U 3 (x 
Ui(x 
U 5 (x 
U e (x 
U 7 (x 
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3. The sequence {W 'Kt(n)}t,n( f or t = 0, 1, 2, • • • ,5) are provided below. 

WKo(n) = C(n,Q) = 1 (empty graph) 
WK^n) = C(n + 1,1) = n + 1 (•) 

WK 2 (n) = C{n + 2,2) = i(n 2 + 3n + 2) (• •) 



WifsO) = — (4n 3 + 18n 2 + 28n+15 + (-l)") 
16 



( Remarks 4 of Section 6 of Bona and JvJ}$) 

WKAn) = — (2n 4 + 12n 3 + 28n 2 + 30n + 13 
16 



-l) n (2n + 3)) 




WK 5 (n) = (12n 5 + 90n 4 + 280n 3 + 450n 2 + 374n + 129 

192 

+ (-l) n (30n 2 + 90?i + 63)) 
2.5 Star Graphs 

The Star Graph St of order t is a tree with t + 1 vertices, one of them of degree 
t and all others of degree 1. 
• • • S& 



a* 



If we let the hub vertex a have value i, then the rest of all the vertices must 
have values in [n — i] . So 



ra+l n+1 t 



ws t (n) = j> + 1 - if = e fc* = E E A (*> i ) c ( fc - 1 + l ' *)] 

i=0 fe=l fc=l i=l 

/ 

=E^( f i*) c '(»+ i + i >*+ 1 ) 



by Fact|n3l3. in Section IO 
Now 



^(1-x) 



1 + 2 



A(t,i) 



A t (x) 

x{i- x y+ 2 ' 



n 



and 



F(s, y) = P( S t)j = g tf = a;(l-x)(l-xe«)' 



t=o 



by Fact IO 4. 



2.6 Cubic Graphs 

The Cubic Graph C is defined as follows: a 




Given four numbers a, b,c,d G [n], e G [n — max(a, b, c)], 
/ G [n — max(a, b, d)],g £ [n — max(a, c, d)], and ft- G [n — max(b, c, d)]. 

Hence 

W^C(n) = (n+ 1 — max(a,b,c))(n + 1 — max(a,b,d)) 

a,6,c,d— 

(n + 1 — max(a, c, d))(n + 1 — max(b, c, d)) 

n n 

= (n + 1 — max(a, b, c))(n + 1 — rnax(a, b, d)) 

d— mai(fl t 6,c)-0 

(n + 1 — max(a, c, d))(n + 1 — max(b, c, d)) 

n d n 

= ^2(J2+ J2 ) (n + l-max(a,b,c)) 

d=0 k=0 k=d+l max(a,b,c)=k 

(n + 1 — max(a, b, d))(n + 1 — max(a, c, d))(n + 1 — max(b, c, d)) 
= EE E (n + l-d) 3 (n + l-fc) 

d— A;— maa;(a,(),c)-fc 
n— 1 n 

+ E E E 

d—0 k—d-\-l max(a,b,c)—k 

(n + 1 — max(a, b))(n + 1 — max(a, c))(n + 1 — max(b, c)) 
:= [/] + [//]• 
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[I] = (n+l)(?i + 2)(n + 3)(157i 5 +138n 4 + 533n 3 + 1074n 2 + 1180n + 560). 
By the Inclusion Exclusion Principle. 

(n + 1 — k)(n + 1 — max(a, b))(n + 1 — max(a, c))(n + 1 — max(b, c)) 

max(a,b,c)=fc 

fc fc 

= 3^(n + 1 - k) 3 (n + 1 - + l) 2 - Z 2 ) - 3^(n + 1 - fc) 4 + (n + 1 - fc) 4 . 

2=0 1=0 

= 7T^™(™ + l)(n + 2)(45n 5 + 357n 4 + 1177n 3 + 1971n 2 + 1638n + 412). 
3360 



WC(n) = [I] + [II] 

= C{n + 8,8) + 26C(ra + 7, 8) + 175C(n + 6,8) + 316C(n + 5, 
+ 175C(n + 4, 8) + 26C(n + 3, 8) + C(n + 2, 8). 

1 + 26x + 175x 2 + 316x 3 + 175s 4 + 26x 5 + x e 

p{c) = a^F ■ 

Let HCd be a d-dimcnsional Hypcrcubic Graph. Then 



P(HC ) 
p{HCx) 
p{HC 2 ) 
P(HC 3 ) 
p(HC d ) 



1 



(1 -xf 
1 

l + 2x + x 2 

1 + 26a; + 175a; 2 + 316x 3 + 175a: 4 + 26a: 5 + x 6 
(1 - xf 



(We conjecture the last one(general case)!) where P 2 d_ 2 {x) is a symmetric 
polynomial of degree 2 d — 2. 

2.7 Octahedral Graphs 

In this section we focus on octahedral graphs. An octahedral graph is the 
Platonic graph with six nodes and 12 edges having the connectivity of the oc- 
tahedron. Let us consider the following graph. 
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•I 

Let OH be an octahedral graph given in the figure above. Let a top vertex 
have value k and a bottom vertex I as in the figure. Then rest of them have 



values in [n — m], where m = max(k,l). 
(n+ 1) X (n + 1) matrix for which b(n, m)ij 
and b(n, m)ij = 1 otherwise. For example, 

"11111" 



S(4,0) 



5(4,3) 



110 

110 









,5(4,1) 



5(4,4) 



Let B(n,m) = (b(n,m)ij) be an 
= if max(i+j, max(i, j) +m) > n 



,5(4,2) = 



1 











1 


1 


1 








1 


1 


1 








1 


1 


1 







































WOH(n) = ^T((m+ l) 2 - m 2 )trace(B(n, m) 4 ) 



m=0 
fe-1 



2 k 



WOH(2k) = ( 2m + l)trace(B(2k, m) 4 ) + ^ ( 2m + 1 )( 2fc + 1 - m Y 



m— 



== (I) + HI) 

k 

WOH(2k + 1) = Y ( 2m + l)trace{B(2k + 1, : 



m=0 



2fe+l 

^ (2m+l)(2fc + 2 



=: (777) + (IV). 



In order to compute the quantities (I) and (III) above, we need to define 
the following. Let D(n) be an (n + 1) X (n + 1) matrix all of whose entries are 
1, and let S(n, m) := (sy) be an (n + 1) x (n + 1) matrix such that s,j = 1 if 
i + j > 2n — m and otherwise. 

For example, 



5(4) -5(4, 2) = 



1 


1 


1 


1 


1" 




"1 


1 


1 


1 


f 
















0" 


1 


1 


1 


1 


1 




1 


1 


1 


1 


1 



















1 


1 


1 


1 


1 




1 


1 


1 


1 


1 



















1 


1 


1 


1 







1 


1 


1 


1 


1 
















1 


1 


1 


1 










1 


1 


1 


1 


1 













1 


1 
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(D(r) - S(r, m)) 4 = (D(rj 2 - D{r)S(r, m) - S(r, m)D(r, m) + S{r, m) 2 ) 2 
=: (M -DS-SD + K) 2 

= M 2 - AIDS - MSD + MK - DSM + DSDS + DSSD - DSK 
- SDAI + SDDS + SDSD - SDK + KM - KDS - KSD + K 2 . 

Let us define 

p(r,m) := trace((D(r — 1) — S(r — l,m)) 4 ) 

= r 4 - 4C(m + 1, 2)r 2 + 4[C*(m + 1,3) + C(m + 2, 3)]r 
- 4C(m + 2, 4) -C(m + 1,2). 

Then we have 

fe-i 

WOH(2k) = J2( 2m + l)p(2fc + 1 - m, 2(k - m)) 

m=0 
2/v 



^(2m + l)(2/c + l-?Ti) 4 

(fc + l)(2fc 2 + 2k + l)(12/c 3 + 30/c 2 + 27k + 10), 



m— k 

l 

~~ 10 

WOH{2k + 1) = ^ (2m + l)p(2fc + 2 - m, 2fc + 1 - 2m) 

m=0 
2fc+l 

+ ^ (2m + l)(2fc + 2-m) 4 

m=k+l 

= — (k + l)(2fc 2 + 6fc + 5)(12/c 3 + 42/c 2 + 51k + 20). 



Finally, we get 
WOH(n) = V" (2m + l)p(n + 1 - m, n - 2m) 

m=0 
n 

+ ^2 (2m.+ l)(n + 1 - m) 4 

= — (6n 6 + 54n 5 + 210?i 4 + 450?i 3 + 559n 2 + 381n + 115 
160 

+ (-l)' l (10n 3 + 45n 2 + 75n + 45)) 

and 

1 + 7x + 48a; 2 + 89a; 3 + 142a; 4 + 89a; 5 + 48x 6 + 7a; 7 + a; 8 



P (OH) 



(1 + a;) 4 (l -x)< 
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2.8 Complete Bipartite Graphs 

Let K Piq = (V 7 E) be a complete bipartite graph of order (p,q). So, V = 
X U Y, E = {ab = ba\a € X, b G Y}, \X\ = p, \Y\ = q. We also let WK p , q , a be a 
weighted graph with a distribution a of the form a = ((jii, n 2l • • • , (mi,m 2 , • ■ 
and r = max(ni, n 2 , • • • , Then rrij G [n — r](j = 1, 2, • • • , g). 



WK p ^(n)= (n + 1 - max(n 1 , n 2 , ■ ■ ■ , %,))'' 

71 1 i n 2 j" ' ,?lp — 

n 

= ^((fc+l) p -F)(n + l-fc)« 
If we let l(t, k) = A(t, k + 1) and A t (x) = Yfk^o ^(*> /c ) xfe - Then 

OO 

p{K v , q ) = \WK p , q (r 

For example, 



A p (x)A q (x) 

v ;l (l-x)P+i+ 1 
n=0 y 1 



A 3 (x) = 1(3, 0) + 1(3, l)x + 1(3, 2)x 2 

= l + ix + x 2 , 
A 4 (x) = 1(4, 0) + 1(4, l)x + 1(4, 2)a; 2 + 1(4, 3)x 3 

= 1 + Ux+ llx 2 +x 3 7 

and 

,„ , / T r \ (l + 4a; + x 2 )(l + lla:+lla: 2 + x 3 ) 
P(K 3A ) = p(K 4>3 ) = (1-s) 8 ' 

Remark 2.9. 

1. We have H(x,y,z) := £~ =o P(Kp,«)§§ = (ra^f^ 
£ FFe ftowe p(K p , p ) = p^pFn • 

5. T/ie graph K p . q has no cycles of odd length. Hence generating functions 
are not quasi-polynomials (as we will show later). 

4- The degree of the numerator in p(K p ^ q ) is p + q — 2 and the numerator is 
symmetric. 

3 Rational Polytopes and Rational Generating 
Functions 

Readers unfamiliar with this topic may wish to consult the book of R. Stanlevf|21|. 
Chapter 4). A quasi (or pseudo)-polynomial of degree d with a quasi-period 
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N is a function / : N — ► C of the form f(n) = J2i=o c i( n ) nt where the coeffi- 
cients are periodic functions of a common period TV and the leading coefficient 
c c i{n) is not identically zero. 



Example 3.1. The number f(n) of unit squares in the region bounded by x 

3' 



0,x = n(> l),y = 0,y = fa; 



|n — 2 if n = Q(mod3) 
f(n) = ^ |n — | if n= l{modZ) . 

\n — \ if ' n = 2(mod3) 



C\{n) — | is constant(l-periodic), but co(n) is 3-periodic. 

The generating function associated with the sequence {f(n)}^L is as follows: 



x 2 (l + x + 2x 2 ) 
(1 -x)(l -x 3 ) ' 



g{x) = f( n >" 



Here is another nontrivial example from the page 220 (Exercises 13.9 and 
13.10) of Pach and Agarwal [T5| . 

Example 3.2. Given a set P of n points in the plane, for any p G P, let pLp{p) 
denote the number of farthest neighbors of p, that is, 

Vp{p) = \{Q 6 P \\P~q\ = max reP \p- r\}\. 

Let fi(n) = max\p\— n X) p ep A*p(p)- It turned out (see Csizmadia 11UI/ or Avis 
et al Wj) that if n is sufficiently large, then 



'.in 



if n= 0(mod2), 



/"(«) = '— + y + \i ifn=l(modA), 



if ifn = 3(mod4). 
The generating function associated with the sequence {/x(n)}5£L * s as follows: 

3 - 3a; + 4a; 2 - x 3 - 3a; 4 + 2a; 5 



9{x) = Y,f J '( 



nix 



(l-x) 2 {l-x 4 ) 



If denominator in the reduced generating function associated with a certain 
sequence {f(n)}^ =0 has a factor 1 — x N then /(n) has a quasi-period N as in 
Example 0(1 - a; 3 ) and Example El (1 - a; 4 ) [2T]. 

Suppose we have a finite set of points {t\, t2, ■ ■ ■ , £d} in M m . The convex 
hull of the set {t\, t-z, ■ • ■ , td} is the set of all convex combinations of the given 
points, i.e. {a; £ R m : x = Yli=i^i> A — 0, i A, = 1}. By a convex 
polytope, or simply a polytope, we mean a set which is the convex hull 
of a non-empty finite set {ti,t2,--- ,td}- An afflne combination of points 
^1)^2) • " i tk from R m is a linear combination Ai<i + \2t2 + • • • + Xktk, where 
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Ai + • • • + Afc = 1 and Aj g M for i = 1, . . . , k. A fc-family (ti, t2, • ■ ■ , ifc) of 
points from K m is said to be affinely independent if a linear combination 
Aiti + \2t2 + ■ • • + Afetfe with Ai + ■ • ■ + Afe = can only have the value when 
Ai = A2 = • ■ ■ = Afe = 0. A polytope V with the property that there exists 
an affinely independent family (ii,<2, ■ ■ ■ , t£) such that V is a convex hull of 
{t\,t2,--- ,td} is called a simplex. (Refer Barvinok Br0ndsted Miller 
and Sturmfels [TJI, Stanley [20j or Schrijver ^U] for details on polytopes.) 
For a graph G = (V, E) g 5, let to = | F| and 

P(G) = {{xx,X2, ■ ■ ■ ,x m ) g |xj +Xj < 1 for all edge g S}. 

Note that for every simple graph G V(G) is a polytope which is contained in 
the m-dimcnsional unit hypcrcube. If all of the coordinates of the vertices of 
the polytope are integers, then we call such a polytope an integer polytope 
(or an integral polytope, a lattice polytope). If all of the coordinates of 
vertices of the polytope are rational numbers, then the associated polytope is 
called a rational polytope. 

Note that every V is homeomorphic to a ball B d , for some d. This d is the 
dimension dim(V) of the polytope V . We denote the boundary (resp. inte- 
rior) of V by ST 5 (resp. V). a £ V is a vertex of V if there exists a closed affinc 
half-space H such that V C\H = {a} 

If V g R m is a rational convex polytope and if n g Z+, then we define 
i(P,n) := \nVnZ m \ and i(T,n) := \nPf]Z m \. 

This is called the Ehrhart quasi-polynomial of V and V respectively. (We 
will show later that this is a quasi-polynomial.) If V is an integer polytope then 
this is simply called the Ehrhart polynomial of V and V respectively. 

If P C Q m , then define den/3 (the denominator of (3) as the least integer 
qeZ + such that qf3 C Z m . 

Theorem 3.3 (Stanley, 1997 |2l], pp237-238). If V is a rational convex 
polytope of dimension d in W 71 , and 

00 

F(V,x) := l + ^i(V,n)x n . 

n=l 

Then F(V,x) is a rational function S^j of x, where deg(P(x)) < e? and Q(x) 

can be written as Y\p^ v {^ — x denl3 ). If F(V,x) is written in lowest terms, then 
x = 1 is a pole of order d+1, and no value of x is a pole of order > d + 1. We 
also have the following (reciprocity for Ehrhart quasi-polynomial) : 

z(7>,n) = (-l) d z(P,-n). 

Theorem l3 . 31 savs that i(P , n) is a quasi-polynomial with the correct value 
i^, 0) — 1, and D{x) = Ili3gy (1 — a; de ™' 3 ) is not in general the least denominator 
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of F(V,x). However, the least denominator has a factor (1 — x) d+1 but not 
(1 — x) d+2 , while D(x) has a factor (1 — See the following example 

(Stanley |3T]). 



Example 3.4. Let V he the convex hull of the vertices set 
V* = {(0, 0, 0), (1, 0, 0), (0, 1, 0), (1, 1, 0), (1/2, 0, 1/2)}. Then 



D(x) = ' [ (1 - x denP ) = (1 - xf{\ + x), 

pev* 



but F(V,x) 



(i-z) 4 



In order to prove our main theorem we need some information about the 
denominator of the vertex set for the polytope V{G). 

Lemma 3.5. Let G = {V, E) with \V\ = m be a simple bipartite graph and 
V* the set of vertices of the polytope V(G). Then den{V*) — {1}. Otherwise 
(equivalently, if the simple graph G is not bi-colorable, or if it has a cycle of odd 
length), then den(V*) = {1,2}. 

Proof. If G = (V, E) is bipartite. Then the defining matrix for the polytope 
V(G) is totally unimodular (Example 1 on p. 273, fH). Thus, den(V*) = {!}. 

If G = (V, E) is not bipartite, then, G contains at least a cycle with odd 
length. Consider a cycle with odd length as a subgraph and a submatrix in the 
defining matrix for the polytope V(G). 

Let A be the defining matrix for G, and suppose the polytope V{G) is defined 
by the system {x £ R m : Ax < b, x > 0} where ylisasxm matrix. Then we 

can rewrite this system as {x £ K m : A'x < b'} where A' = ^ ^ ^ where 

L m is an identity matrix of order m and b' = ( ^ J where is a m dimensional 

vector with all zeros. 

Then note that any vertex of the polytope V(G) is determined by the unique 
solution of a subsystem 

A"x = b" 

where A" is a m x to minor of A' with &et(A") ^ and 6" is a to sub- vector of 
b' (Theorem 8.4 Q]5]). Thus all we have to show is that det(A") = {±1, ±2} if 
&et(A") ^ 0. Then, there are three cases we have to consider. 

1 . All rows of A" are from I m . 

2. All rows of A" are from A. 

(a) A" contains a cycle (or cycles) with odd length. 

(b) otherwise. 

3. k rows of A" are from A and m — k rows are from /„, . 
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(a) A" contains a cycle (or cycles) with smaller odd length k. 

(b) otherwise. 

If A" = I m , then clearly this defines the origin. For Case 12 bl and Case 13 bl A" 
does not contain a cycle with odd length. Thus we are done. For Case I2al and 
Case |3al we have to show that the matrix for each subsystem of the system 
A" x = b" defining a cycle with its length k, where k is an odd positive integer 
with k > 3, has its determinant ±2. Let A be the defining matrix for a cycle 
with odd length. Note that A is a k x k matrix and note that after permuting 
columns and rows, A forms such that: 



/ 1 1 
1 1 




\ 1 



\ 




1 1 
1/ 



Note that A is an upper triangular matrix except for the 1 in the first position 
of the fcth row. So the determinant of A is clearly 2. Thus, den(V*) = {1,2}. □ 



The next theorem is the conclusion of all that we discussed. Its proof follows 
immediately from Theorem 13.31 and Lemma f3. 51 

Theorem 3.6. If a simple graph G = (V,E) with \V\ = m is bipartite (equiva- 
lently, x(G) then WG(n) = i(P(G),n) is an Ehrhart polynomial and 



P(G) 



P(x) 



(1 - x) m + x ' 



where P{x) is a symmetric polynomial of degree < m. That is, ifWG(n) is an 
Ehrhart quasi-polynomial in n of a quasi-period ^ 1, then the quasi-period is 2, 
G contains a cycle of odd length (hence x{G) > 2J and 



P (G) 



P(x) 



(l-x 2 ) k (l-x) m + 1 - k ' 
where P{x) is a polynomial of degree < m + k and k is a nonnegative integer. 



Corollary 3.7. If a simple graph G = {V, E) with \V\ — m is either a tree, 
a circular graph of even length, a discrete graph, a hypercubic graph, a com- 
plete bipartite graph, or a grid graph, then p{G) = ; where p(x) is a 
symmetric monic polynomial of degree at most m. 

Proof. All of those graphs are bipartite graphs. Hence, by Theorem 13.61 the 
result follows. □ 
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Example 3.8. Consider the Circular Graph C 3 (a cycle of odd length 3). 
P(<3$) = {(r,s,t) G R 3 |r,s,t > 0, r + s, s + t,t + r < 1} 

v* = {(o, 0, o), (1, 0, o), (o, 1, o), (o, 0, 1), (i i i)} 

m = dim(V(C 3 )) = 3 

^C 3 (n) = i(P(C 3 ), n) = T-(4n 3 + 18n 2 + 28n + 15 + (-1)") 

Id 

,«?,) = F {nCM = l + ±mC 3 ),n^ = (1 ^ 2 ) ( t! 2 . )3 = (I^^i 
= JJ (1 x den/3 ) = (1 - x- 2 )(l - x) 4 . 



Example 3.9. Consider the Circular Graph C4 (a cycle of even length 4)- 
ViCi) = {(r,s,t,u) G M 4 |r,s,i,u > 0, r + s, s + t,t + u, u + r < 1}, 

y* = {(0,0, 0,0), (1,0, 0,0), (0,1,0,0), (0,0, 1,0), (0,0, 0,1), (0,1,0,1), (1,0, 1,0)}, 

m = dim{V{Ci)) = 4, 

WC 4 (n) = C(n + 2, 4) + 2C(n + 3, 4) + c(n + 4, 4) = 1 + + \n 2 + n 3 + -U 4 , 



P (C i ) = l + J2K'P(C i ) 



n)x 

n=l 



l + 2x + x 2 P(x) 



(1-x) 5 {1-XY+ 1 ' 



D(x)= ' I (1 - x denl3 ) = (1 - x) 7 . 

f3<£V 



Example 3.10. Consider the simple graph G — (V, E), where V = {1, 2, 3, 4, 5} 
and E = {12,23,34,45,15,13}. Hence, the associated polytope and its vertices 
are (we used LattE Mty and CDD \14Jj for computational experimentation): 

V(G) = {(r,s,t,u,v) G K 5 |r, s,t,u,v > 0,r + s,s + t,t + u,u + v,v + r,r + t < 1}. 

V* = {(0, 0, 0, 0, 0), (1, 0, 0, 0, 0), (0, 1, 0, 0, 0), (0, 0, 1, 0, 0), (0, 0, 0, 1, 0), 

(0,0, 0,0,1), (1,0,0, 1,0), (0,1, 0,1,0), (0,1, 0,0,1), (0,0, 1,0,1), 

(1/2, 1/2, 1/2,0, 0), (1/2, 1/2, 1/2, 1/2, 0), (1/2, 1/2, 1/2,0, 1/2), 

(1/2,1/2,1/2,1/2,1/2)}, 

dim(T(G)) = 5, 

„,_, . 121 535 219 . 13 , 45 , 3 5 . .... 7 3 1 2 , 

WG(n) = 1 n H n 2 H n 3 H n 4 H n 5 + f-l)"( 1 n H n 2 ), 

v ; 128 192 64 6 64 32 v ' 428 64 64 h 

1 + 7x + 22.T 2 + 30a; 3 + 22x 4 + 7x 5 + x 6 

P( - ' ~ (l-x) 3 (l~x 2 ) 3 ' 

D{x) = Y[ (1 - x den0 ) = (1 - x) w {l - x 2 ) 4 . 
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C3 in Example 13.81 has a cycle of odd length (3) and C4 in Example 13. 91 has 
a cycle of even length (4). However, the graph G in Example 13.101 has three 
cycles of length 3, 4 and 5, two of them odd. 

Remark 3.11. If a simple graph G = (V,E) with \V\ = m is bipartite and 
p(G) = ^ P ^l+i , then coefficient vector h(G) = (ho, hi, ■ ■ ■ ,h m ) of a symmetric 
polynomial P(x) = ho + hi satisfies the following properties (see 

Stanley 

(1) h = 1. 

(2) h m = (-1)"WG(-1). 

(3) min{j > 0\WG(-1) = WG(-2) = ■■■ = WG(-(m - j)) = 0}. 
= max{i\hi 7^ 0} 

For example, WC 4 (-1) = WC 4 (-2) = and h(C 4 ) = (1,2,1,0,0). 

(4) WG(—n — k) = (— l) m WG(n) for all n if and only if hi = h m +i-k-i for all 
i and h, n+2 -k-i = h m +s-k-i = ■ ■ ■ = h m = 0. 

(5) hi > for all i (Nonnegativity) (See Stanley XEkXHSj .) 

(6) If Gi is a subgraph of G 2 , then V{G X ) C V{G 2 ) 

and hi(Gi) < hi(G 2 ). (Monotonicity)(See Stanley JMffl!/ . ) 

(7) volumeifPiG)) = — leading coefficient ofWG(n). 
(See Stanley \T3\W&j .) 

(8) All real roots a of WG{n) = satisfy -m < a < |_f J . (See Beck et al 

Remark 3.12. There seem to be some connections to semi magic cubes. A 

semi magic cube is an k x k x • • • x k table with nonnegative integral entries 
such that each directional sum must be equal to n (a magic sumj. (The reader 
can find more details on semi magic squares and semi magic cubes in JhT/). For 
example, suppose we have a 2 x 2 magic square. Then, the defining polytope for 
a 2 x 2 magic square with a magic sum n is a face of the polytope defining a 
weighted cycle graph WCA(n). In general suppose we have a 2 d magic square. 
Then the defining polytope for a 2 d magic square with a magic sum n is a face 
of the polytope defining a weighted d— dimensional Hypercubic Graph HCd- 

Remark 3.13. We can get p(G) from the information of the graph G = (V, E) 
using the package LattE by J.A.De Loera et. al. (JTTj ), and the Elliott Maple 
package by G. Xin (I24!j ), which improved the Omega package by G. Andrews 
et. al. (Jl]l or references therein (in fact, their serial articles)). We also can find 
the coordinates of all vertices of the polytope P{G) using edd and cdd+ by K. 
Fukuda (\T^). 

4 Further Questions 

1. It is obvious that, for a given number m of vertices, 

WK m {n) < WG(n) < WD m (n) 

for any simple graph G = (V,E) with \V\ = to, where K n (resp. D m ) 
is a complete(rcsp. discrete ) graph of order to. We can ask the same 



22 



question for trees. That is, given the number of vertices what kind of 
trees T achieve the maximal or minimal values of WT{n). 

2. If a simple graph G is connected and p(G) = ^fj, where P(x) and Q(x) 
are polynomials of lowest degree, then is P{x) always symmetric ? That 
is, if k = degree(P(x)) then does x k P(^) = P{x) hold? 

3. For a given simple graph G, can we express the volume of the polytopc 
V{G) in terms of the graph G? 

4. If • is an operation between two simple graphs Gi, G2 and it is closed in 
the simple graphs, what is W{G\»G2){n) and p(G\*G2)7 Can we get any 
relation between WG(n) (resp. p{G)) and WG(n) (rcsp. p(Gj)? (G is a 
complement graph of G). 

5. Given two simple connected graphs G\ = (V,Ei) and G2 = (W, £2), 
choose a vertex v <G V,w <G W . Make v and w adjacent by adding an 
edge between them so that two separate graphs G± and G2 becomes one 
connected graph G. Which vertices in each side do we have to choose in 
order to maximize or minimize WG(n) ? 

6. Generalization or Extension of the Weighted Graph by inserting one or 
several slack variables in the edge between two adjacent vertices(cf. Gear 
Graph), or to the Problems related to (or symmetric) Magic Squares. 

i • • 3 ^ 



AI = (fain, k)) : (n + 1) x (n + 1), where 

fa ( n > fc ) = an-.-,- c ( n + k-i-j,k) 

for < i,j < n. Note that M = (/y(n,0)) = B(n). 

7. For what kinds of rational function f(x) does there exist a simple graph 
G so that p(G) = f(x)? In other words, what is the image p(Q) of the 
map p in Z[[x]] ? 

8. Compute WG(n) and p(G) for other shape of simple graphs, like Wheel 
Graph, Cayley Graph, Complete k-partitc Graph, Web Graph or Grid 
Graph, Regular Graphs, other Platonic Graphs etc... 
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